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We extend the application of axial Ward identities to calculate bA, bp and br, coefficients that give the mass dependence of the 
renormalization constants of the corresponding bilinear operators in the quenched theory. The extension relies on using operators 
with non-degenerate quark masses. It allows a complete determination of the 0{a) improvement coefficients for bilinears in the 
quenched approximation using Ward Identities alone. Only the scale dependent normalization constants Z% (or Zg) and Zt are 
undetermined. We present results of a pilot numerical study using hadronic correlators. 
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Wilson's discretization of the Dirac operator introduces 
lattice artifacts at 0{a), the effects of which are large for 
lattice spacings accessible in present simulations. It is thus 
expedient to devise improved lattice discretizations, and 
significant recent progress has been made in this direction. 
In particular the ALPHA collaboration has implemented 
Symanzik's improvement program in which both the 
action and external sources are improved by the addi- 
tion of higher dimensional operators. A key ingredient is 
the development of methods to determine the coefficients 
of these extra operators, the "improvement coefficients", 
non-perturbatively. In this note we present an extension 
of these methods which allows the determination of all 
improvement coefficients for bilinear operators. 

We begin by reviewing the results of the ALPHA col- 
laboration. They have shown that 0(a) artifacts can be 
removed from on-shell quantities by the addition of a sin- 
gle local operator of dimension five resulting in the 
Sheikholeslami-Wohlert (or "clover") action p] 
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Ssw = a^csw ^'ip{x)j<Jf,^Ff,^{xyj{x) 



(1) 
(2) 



Here Sg and Sw are respectively Wilson's gluon and 
quark actions, and Ssw is the Sheikholeslami-Wohlert 
term {F^^, is the lattice gluon field strength tensor). Im- 
provement of quark bilinears is slightly more involved . 
We consider only flavor off-diagonal operators, the bare 



lattice forms of which are 

-0-' 7^75-0'' , 



4(ife) 



^-^7^0 , 



(3) 
(4) 
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with j ^ k being flavor indices, cr^i. = [7^,7y]/2, 75 = 
71727374, and the hermitian 7 matrices satisfy {7^, 7^} = 
2(5^1^. Removal of 0{a) errors from the on-shell matrix 
elements of these bilinears requires both the addition of 
extra operators (except for P and S), 



A, 



and the introduction of mass dependence, 



(8) 
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Here X ~ A,V, P, S, T, the are renormalization con- 
stants in the chiral limit, and rriij = (mi + mj)/2 is the 
average bare quark mass. 

Complete 0{a) improvement of matrix elements re- 
quires that the coefficients cswi cx, and bx, as well 
as the matching constants Z^, be determined non- 
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perturbatively.0 Previous work has shown how the en- 
forcement of axial and vector Ward identities (WI) allows 
one to determine Z^, and Zp/Zg csw, CA and by 
cv i, CT and bp - bA and bs We dis- 

cuss here an extension that yields bA, bp, and &t-0 This 
provides an alternative to the non-perturbative method 
proposed in Ref. which uses the short-distance behav- 
ior of two point functions. The two remaining constants 
Zp (or Zg) and Z^ are scale and scheme dependent, and 
cannot be determined using WI. 

It is important to note that the relations we derive do 
not extend directly to the unquenched theory, which re- 
quires additional improvement constants and a more com- 
plicated set of conditions as will be presented in |16|. 

We begin by recalling the ALPHA method for determin- 
ing Csw and ca The improved axial current should 
satisiy 



(12) 



when inserted between on-shell states. Here is the 
renormalized quark mass. It follows that the ratio 



2to, 



(mf' + mf ) 



p Z^(l + bparriij) 



^ ' Z\{\ + bAam,,) 



(13) 
(14) 



should be independent both of the choice of sources J and 
of the time t, as long as t ^ 0, up to corrections of 0{a?). 
This is to be achieved by simultaneously tuning csw and 
Ca ■ Our implementation of this condition differs from the 
Schrodinger functional method of Ref. in that we use 
standard two-point correlation functions, with a variety of 
choices of sources. We also fix csw « priori and use the 
condition to determine only ca- 

It is convenient for the following discussion to introduce 
b coefficients defined using the masses fhij , which we refer 
to as WI masses, rather than the bare quark mass, i.e. 



(Xfl)(»^) = Zl{l + bxarh,,){Xif^^ 



(15) 



The advantage of the WI mass is that it can be deter- 
mined, using eq. ([l3|), without the need for chiral extrap- 
olation. By contrast, to determine the bare quark mass, 
am = 1/(2k) — 1/(2kc), one needs to know the critical 
hopping parameter, Kc, the calculation of which requires 
chiral extrapolation. The difference is particularly signifi- 
cant when using standard two-point functions as opposed 
to the Schrodinger functional. To determine the standard 
bx from the bx we need the relation between the bare 

^For brevity we refer to the as improvement coefficients in the 
following. 

■^A preliminary account of this work was given in Ref. pa]. 



and WI masses, which is given in eq. 
quenched theory it follows that 



Z° 

= -^oj^bx 

-^A^S 



0{a) 



below. In the 



(16) 



Note that we only need this relation at leading order in a, 
since the bx appear only in 0{a) corrections. 

With Csw fixed, Zy and by can be obtained in the stan- 
dard way using charge conservation. We use the forward 
matrix elements of (V/)4 between pseudoscalars, 
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Z^(l + 6y 07712) 



(^.P(i2)(f,r)J(2i)(0)) 



(17) 



with T > t > and J = P oi A4. Here and below we 
use an abbreviation for the WI mass for two degenerate 
flavors: rhi = rha. Note that the cy term in V/ does not 
contribute to the r.h.s., and so is not determined by this 
condition. 

The remaining improvement constants are determined 
by enforcing the generic form of the integrated axial WI 
(AWI), in which the operator O^^ = VS^^^PV'^^) is trans- 
formed into 50^^ = ^^^^^^Ttp^^") by a chiral rotation in 
the 1, 2 flavor subspace: 



Here SS is the variation in the action 



(mf-Kmf)(P«)(i2)-a,(AH)(i^) 



, (19) 



and the chiral rotation is restricted to a 4-volume V con- 
taining y but not z. The Ward identities should hold up 
to corrections of O(a^) if the operators and action are ap- 
propriately improved. There is, however, an obstacle to 
implementing these constraints, arising from the integral 
over X in SS. This brings the pseudoscalar density Pp{x) 
into contact with the operator Op{y) on the l.h.s. of (p^), 
implying that on-shell improvement of P and O is insuffi- 
cient to improve the AWI. (The integral over d^A^ gives 
a surface term which does not involve contact with O.) 
The problematic term is explicitly proportional to quark 
masses, and thus is absent in the chiral limit. For this rea- 
son the AWI has previously been used only in the chiral 
limit, from which one can determine Z^, Zp/Zg, cv and 
ct, but not the bx- 

Our new observation is that by looking at the detailed 
dependence on the quark masses one can work away from 
the chiral limit and yet avoid contact terms. This allows 
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the determination of certain linear combinations of the 
bx- The cost is the need to use non-degenerate quarks. 

For all but the coefficient br, it is sufficient to work 
with two non-degenerate quarks: toi — m2 7^ m^. Since 
the contact term is explicitly proportional to + m|^, 
which itself is proportional to rhi2 at small quark masses, 
one can remove it by extrapolating to rhi — rfi2 ~ 0. 
The remaining 7713 dependence allows one to determine 
combinations of the bx- In the following, we describe 
how this works for the different bilinears, and explain our 
particular implementation. The extrapolation to rhi = 
7712 = will be implicit throughout this discussion."^ 

As a first application we show how to obtain bA — by, 
as well as cy, using the AWI with O = Vf^. The identity 
for the time component can be written 



ri 



12: 



E^((^/)f'(y,y4) J(3i)(o)) 

^0(1 + bAam3/2) 



bva,mj/2) 



(20) 



where 5Si = — /y df^{Ai)^. For the source we take J = P 
or A4. We have chosen to put the current at p = so 
that the cy term in Vj does not contribute. Since we 
know csw and ca we can calculate ri , and determine bA — 
by from its dependence on 7713. The intercept provides a 



determination of Zy independent of that from eq. (17) 



The AWI for the spatial components can be written as 



Y.y.e-p-y{5sf^^ [x^. + acya^r,^](23)(y,y4) (0)) 



(31) 



(21) 



Enforcing this equality (for any small ih^) provides a de- 
termination of Cy. 

Reversing the roles of vector and axial bilinears provides 
a second determination of bA — by. For example, if cy is 
known from above, then one can use the ratio 



Ege'^-g(j5f^^ {Ai)r'{y,yA) vrm 

Z°{1 + byam3/2) 



(23), 



.(31) 



r2 



(22) 



Z% ■ Z\{1 + bAam3/2) 

Given Zy, this also yields Z^. The same information can 
be obtained from the combinations 
1 



7" 



(23) 



^ In practice we keep the term J {^i +"^2)(^'ij)^^^' in eq. ( |l8[ ) prior 
to extrapolation, since this increases the range of to for which the 
ratio of correlation functions is constant. 



l- = ^{l + ibA-by)arh,/2). 

^2 ^y 



(24) 



Applying the method to the AWI with O = P and S 
gives two determinations of bp — bs- For example, one can 
use the rh^ dependence of the ratio 
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J:.e^p-y{pm{y,y4) J(3i)(0)) 
Z^{l + bparh3/2) 



Z°^- Z°g{l + bsam3/2) 



(25) 



with J = F or A4. Given Z^, this also yields Z'^/Z^ 
Alternatively, one can use 
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r4 



j:^^e^P-y{SS\''^ P^^^\y,y4) J^^'HO)) 

S.e^?-?(5(13)(y,y4) j(31) (Q)) 

Z°{1 + bsam3/2) 



• Z0(H-6pam3/2) 



(26) 



where for the source we choose either J*^"^^-* — S''^^-' or 

J2^P^^^^Hz, Zi)P^^^HO) for < y4 < Zi- 

The method described so far does not work for the ten- 
sor bilinear, because the chiral rotation transforms it back 
into (other components of) itself, and the dependence on 
bx cancels. One can, however, use the method to deter- 
mine CT |l3| . For example, the AWI with O 
p = and J — Tk4, can be rearranged into 



Tij at 



1 + acT 



Ey{Ttl'\y,y,)Tl:f\0)) 

EyiSSi''^ Tl''\y,y,)Tl:l'\Q)) 
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(27) 



Here we have moved the ct dependence in (r/)fc4 onto the 
l.h.s., and used the fact that (Ti)ij has no contribution 
from the ct term at p = 0. Given 2'^, this equation de- 
termines Ct- a consistency check is that the result should 
be independent of 7713 . 

It turns out that one can determine bx using the AWI 
(p^), but to do so one must work with non-zero mi and 
7712. This means that the AWI is not completely improved: 
terms of 0{a) result from the contact of the pseudoscalar 
density in SS with the operator O. The key point, how- 
ever, is that these terms are proportional to 77ii + 7712, 
while the dependence on br is proportional to the differ- 
ence TTii — 7712- By separating these two dependences one 
can, in principle, determine bx- 

To explain this in detail we recall that off-shell 0(a) 
improvement requires the addition of an extra operator 
(multiplied by an extra improvement coefficient) for each 
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bilinear [y_3| . For the pseudoscalar and tensor the required 
additions are 



p(12) , p(12) 



+ a4V'^^^75(^ + m2)V'(2^ (28) 
+ ac'rpi^^^haf.^il/) + m3)7/'(3) .(29) 



Inserting the ofF-shell improved operators into the AWI 
considered previously we find 



j:.e^P-y{SS\''^ {Ti)f\y,y,) T^f ^(0)) 

E,-e'^^-^-((r,)if (y,2;4)rif)(0)) 
1 + abrirhi — 1112)/ 2 



Z\{1 - abAmi2) 



— 2a{c'rp + c'p)rhi2 



(30) 



Note that both the bA term and the contact terms pro- 
portional to c'rp + c'p depend on mi2. Thus, assuming that 
ct is known, br can be obtained from the dependence of 
rs on mi — 1712- Note that the ratio depends on ct at 
both zero and non-zero momentum. 

A similar extension can be considered for the other 
AWI discussed above. It is straightforward to see that 
the rhi — TO2 dependence of the ratios ri and r2 allows a 
determination of by + bA, while that of and gives 
a determination of bp -\- bs- Thus, in principle, one can 
determine all the bx using the generic AWI. One can also 
determine the five additional improvement constants c'-^ 
using the dependence on as will be discussed in [ p^ . 
Further consistency checks are provided by the three-point 
vector WI with non-degenerate masses — although these by 
themselves do not allow one to disentangle the bx from 
the c'x- 

Petronzio and di Divitiis have shown that one can also 
use the two-point versions of the vector and axial WI to 
determine a subset of the quenched bx, namely bA — bp, 
bs and by |14|. The key point is again the use of non- 
degenerate quarks. In our numerical study we use some 
of their results, which we recall here. 

The first result is obtained by comparing the WI mass rh 
to the bare quark mass m. These two masses are both re- 
lated to the renormalized quark mass, the former through 
eq. (0), and the latter by 



mf = Z^(l -I- bmarhi) rrii 



(31) 



Combining these relations, and considering only degener- 
ate masses m; = mi, one finds 



m 
m 



bm)arh\ 



(32) 



^^[l + (bA-bp + bs/2)am] 



(33) 



To obtain the second line we have used the results Zll^Zg = 
1 and bs = —26m, the latter valid only in quenched QCD 
fl^ . Thus, from the slope and intercept of (33) one can 
determine bA~bp + bs/2 and Z'^Zg/ Zp, respectively. To 
use this method we need to determine Kc, which involves 
a chiral extrapolation. 

The relationship between m,;j and {mi + mj)/2 for non- 
degenerate quarks gives additional information. It turns 
out that this information can be gleaned without reference 
to the bare quark masses, and thus without the need for 
chiral extrapolation. In particular, by enforcing (2toi)^-|- 
(2m2)^ = 2(toi -I- 1712)^ to 0(a), one finds 



bA 



4mi2 - 2[mii -|- TO22] 
a[mii - 77122]^ 



(34) 



This result, in terms of bare masses, was already noted in 
Ref. 0. 

The final result from Ref. uses the vector two-point 
WI. Requiring ^^(Vfl,)^,^^^ = (mf - m§)s'-p'^^ leads to 



^12 



^^e^p.^(^(l2)(2?^^)J(2l)(o)) 

^^Ii±l™i(mf-mf). 
Zv[l + 6yami2] 



(35) 
(36) 



We implement this using two sources: j'^^^ — S^^^\ and 
J(2i) = '£sP^'^^K^,Z4)P^^^H0) with <t < Zi. Enforc- 
ing the relation (2mi)^— (2to2)^ = 2(7711 — 7712)^ to 0(a), 
with the l.h.s. determined from the AWI ( p^ and the rhs 
from the vector WI (|35|), we find 



^^^^ 



{bp - bA) = 



A12 - -Rz[777ll - 77122] 
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(37) 
(38) 



We can use eqs. ( |3^ , |3^ ) to determine bp and 65 separately, 
since by and bA are already known from eqs. (|T^,EQ). 



We have performed a pilot test of our method on an en- 
semble of 16^x48 quenched lattices at P — = 6.0. We 
use the tree-level tadpole-improved value for the clover co- 
efficient, csw = 1.4755, rather than the non-perturbative 
value csw = 1-769 ||]. This imphes that our results for 
the bx differ from the non-perturbative values by correc- 
tions of 0{g^). We do not have data with three non- 
degenerate quarks, and so can test only the simpler ver- 
sion of our method which does not yield bx- Nevertheless, 
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our results should suffice to assess the practicality of us- 
ing WI with non-degenerate quarks to determine the bx 
non-perturbatively. 

Previous determinations of the improvement coefficients 
cx have used Schrodinger functional boundary conditions 
in the time direction, with sources J placed on the bound- 
aries. One advantage of this approach is that one can 
work directly in the chiral limit. In our study we use the 
same correlation functions as used in studying the spec- 
trum and decay constants, i.e. we have periodic boundary 
conditions in the time direction. Thus a secondary output 
of our study is a comparison of these two approaches for 
determining the cx- We stress, however, that our method 
for determining the bx works for any choice of sources J, 
and in particular can be applied to the Schrodinger func- 
tional. 

The correlation functions required in the integrated ax- 
ial Ward identity are obtained as follows. Quark 
propagators are calculated using a Wuppertal smeared 
source at t = for five different values of n corresponding 
to aM^ = 0.57,0.50,0.43,0.35,0.24. These propagators 
are used both to construct two-point functions and also as 
sources for propagators with the insertion of SSi defined 
in eq. (|9|). Our insertion volume V is the region between 
t = 4 and 18. The second inversion uses the same k, 
so that, as already noted, mi = m2. To construct three- 
point functions, propagators with and without sources are 
contracted to form the operator O. This allows us to in- 
sert any momentum into O and to place it anywhere in 
the interval 4 < i < 18. 

For the chiral extrapolations we ignore the correlations 
in the data between the different mass points. We find 
poor signals in some of the correlators containing quarks 
of the lightest mass, and so exclude the lightest mass from 
the extrapolations. The remaining four values of quark 
mass correspond to the range rris — 2.6ms, where rus is 
the physical strange quark mass. Note that the relevant 
expansion parameter here is m^a, and this is small for our 
range of quark masses. 

Our results from the various WI are summarized in Ta- 
ble |l|. Where there is a choice of sources, J, we have 
quoted the results with the best signal. In the following 
we discuss the various determinations pointing out salient 
features. 

We determine ca by requiring that the right hand side of 
eq. (13) is as close to a constant as possible over a range 



of times imin — imax- For long times only the pion con- 
tributes and the ratio is constant for all values of ca — thus 
any choice of tmax in the pion-dominated region is equally 
good. To maximise our sensitivity to ca, we choose ^min as 
small as possible while avoiding contact between SSj and 
J. We find that our results are insensitive to small varia- 



eq. 


observable 


two-point d 


three-point d 


E 


CA 


-0.005(15) 


-0.025(10) 


E 




-h0.745(l) 


+0.746(1) 


E 


by 


-1-1.57(2) 


+ 1.57(2) 


E 




-1-0.753(6) 


+0.759(6) 


E 


bA - bv 


-0.56(9) 


-0.48(10) 


Hi 


cv 


-0.66(27) 


-0.65(25) 


Hi 




-M.41(12) 


+ 1.38(11) 


Hi 


bA - by 


+1.2(1.2) 


+ 1.8(1.0) 


||ii 


Cv 


-0.30(26) 


-0.31(24) 


E 


Z'a 


-1-0.76(2) 


+0.77(2) 


H 




-1-0.97(5) 


+0.98(4) 


H 


bA-bv 


-^0.38(74) 


+ 1.09(69) 


H 




+0.95(2) 


+0.98(1) 


Hi 


bp~bs 


-0.09(9) 


-0.06(9) 


1 


CT 


+0.17(7) 


+0.18(6) 






+0.97(1) 


+0.97(1) 


1 


bA-bp + bs/2 


+0.48(1) 


+0.50(1) 




bp-bA 


-0.20(9) 


-0.13(8) 


E 


bs-bv-2bp-2bA 


+0.04(48) 


-0.27(38) 



Table 1 

Results for improvement coefficients from the listed WI. 
The two columns of results are for two choices of dis- 
cretization of derivatives in eq. (p^, as discussed in the 
text. The labels "i"-"v" are explained in the text. 



tions in tmin. There is a similar insensitivity to the choice 
of source J. This is in marked contrast to results from 
the unimproved action [csw — 0) fl^ and gives us confi- 
dence that our tadpole-improved action has small enough 
0{g^a) errors that we can carry out our tests. 

To estimate the size of 0{a^) errors (which become 0{a) 
errors in ca), we use two choices of discrete derivatives in 
eq. ([1^): a scheme based on two-point derivatives, 



dfix + 1/2) = 
d^f{x + l/2) = 

f{x + l/2) = 



f{x + l)-fix), 

f{x + 2) - fjx+l) - f{x) + f{x - 1) 
2 

f{x + l) + f{x) 



and one based on three-point derivatives, 
df{x) - {.f{x + l)- f{x-l))/2, 
d'fix) = f{x + l)~2f{x) + f{x~l). 
The O(a^) correction in the three-point df is four times 
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larger than in the two-point case. The results for ca from 
these two discretization schemes differ at 0(a).^ As shown 
in Table 1^, this difference is substantial as, unfortunately, 
are the statistical errors. Because of this difference, we 
present results for the remaining improvement constants 
using both choices of ca- 

Both discretization schemes yield results for ca signif- 
icantly different from the value ca = —0.083(7) obtained 
previously at csw = 1-769 using the Schrodinger func- 
tional This difference is presumably a combination 
of 0{g^) and 0{a) effects. We are unable to resolve this 
discrepancy in this work. 

The signal in the vector Wl ( |l7| ) is very good, and pro- 
vides our best estimates of Zy and by- There is a tiny 
dependence on ca due to the determination of the quark 
mass using eq. 

We now turn to results from our new method. We use a 
two-point discretization of the derivative in SSi through- 
out.! ^ig- 

we show our results for the quantity ri 
which appears in the AWI (pO[). The data show a linear 
dependence on 7713 for our range of masses. 

There is a numerical subtlety in the extraction of bA — 
bv- When using eq. |2G, we have to make two choices 
concerning the form of Ay^^ : whether to discretize it using 
two- or three-point derivatives; and whether to use the m 
dependent value of ca obtained from eq. |l^, or the chirally 
extrapolated value. Both choices only effect the result for 
bA — bv at 0(a). In particular, it is straightforward to 
see that the two options for ca lead to results for bA — by 
differing by ~ O'M^ jra (assuming pion domination of the 
correlators). Since M'^/m is a large scale, perhaps as large 
as 5 GeV, these differences, although technically of higher 
order, can be numerically significant. We find that they 
are only a 15% effect for the two-point discretization of A/, 
but are much larger for the three-point discretization. For 
this reason we use the two-point discretization. As for the 
choice of c^, we take the chirally extrapolated value, since 
this is the consistent choice at the order of improvement we 
are working. This does, however, have the disadvantage 
of giving poorer plateaus in the ratio ri 

The results for Zy and bA — by are given in Table |. 



^Since our action is only perturbatively improved, one might worry 
that the two determinations of could differ by 0{g^) in addition 
to 0{a). This is not the case, however, because we use essentially the 
same matrix elements in both determinations, and so the difference 
between them is explicitly proportional to a. 

^The results for improvement constants are expected to be fairly 
insensitive to this choice of discretization because, after integration 
over V, different choices differ only by surface terms. 
®We have a similar choice for the appearing in ^sj^^' , and here 
we use the mass dependent value which leads to better signals. This 
does not, however, directly affect our results for the bx since we 
extrapolate to rhi2 = 0. 
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Figure 1. The ratio ri as a function of 1713/2, after linear 
extrapolation to rhi — rh2 — 0. The intercept and slope 
give l/Zy and Sa — by, respectively. The fit excludes the 
lightest point as the expected plateau in the fits to the 
different are not uniformly good. 

Zy is consistent, at the two a level, with the result from 
the vector WI, but has much larger errors. bA — by is 
determined with a 20% error. 

The determination of cy is more problematic and we 
have attempted a number of approaches. The two best are 

(i) We obtain cy by first equating the right hand sides 
of Eqs. ( [20|j2f| ) for each 7713 and then extrapolating 
to TO3 = 0; 

(ii) For each rhi ~ 7712 and 7713, we solve rir2 ~ r^r^, 
where only r2 depends on cy. The chiral extrapo- 
lation in 77ii = 777-2 removes the contribution of the 
contact terms. Lastly we extrapolate to 7773 = 0. 

All methods require knowledge of ca- The correction 
term proportional to cy is ~ 4% in eq. ( pl| ) and ~ 20% in 
eq. (^ , so one needs high statistics to get an accurate es- 
timate. We prefer (i) as it is the most direct, and because 
74 does not have a good plateau. The large uncertainty 
in Cy extracted by these methods accounts for ~ 50% of 
the errors in Z'^, and subsequently in Z% / Z% and ct- 

With Cy (and Zy') in hand, we use eq. ( P3| ) to determine 
our best estimate of Z\. The results from eqs. (p^ ) and 
( [2^ ) are consistent but have larger errors. The estimates 
of bA — by from these two equations have very large errors. 
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Figure 2. The ratio ra as a function of m3/2, after lin- 
ear extrapolation to rhi = rfi2 = 0. The intercept and 
slope give Zp/{Z^Z\) and bp — bs, respectively. The fit 
excludes the lightest point. 



The final application of our new method is the deter- 
mination oibp — bs using eqs. (|25| ) and (p6|). Note that 
neither of these requires knowledge of cy. We find a good 
signal in the ratio but not in r^. The former is shown in 
Fig. H and the slope yields the value quoted in the table. 
The intercept gives an estimate of Zp/ZgZ'^. 

The last mixing coefficient ct is determined from 
eq. (p7|). There is a good signal in all correlation func- 
tions and only ca and Z^ are needed beforehand. The 
error is dominated by the uncertainty in cy which feeds 
m through Z°^. 

Further information, and consistency checks on the pre- 
vious results, are provided by the WI using two-point func- 
tions. Equation ( p3| ) gives a result for Zp/ZgZ'^, consis- 
tent with that from above, with similar errors. It also gives 
a very accurate result for the combination bA — bp + bs/2, 
with little dependence on ca- There is, however, an ad- 
ditional uncertainty due to the choice of Kc. The two re- 
maining equations, ( ^4|) and (^t]), which do not require Kc, 
give rather poor determinations. An important technical 
point is that the 0{a) correction in eq. (^) is large when 
the source J^^^' is chosen to be S'^^^' as MscaiarO. « 1. 
In order to get results consistent with those from the 
choice J(2i) = Y.zP^'^^Hz,Zi)P^^^H'^), we needed to use 





LAiNL 


A T T~)TU A 

ALrMA 


Pert. in. 


Csw 


1.4755 


1.769 


+ 1 




Zy 


+0.745(1)(1) 


0.7809(6) 


-1-0.810 






+0.76(2)(1) 


0.7906(94) 


-1-0.829 




Zjp/Zjg 


+0.77(4)(1) 


N.A. 


-1-0.949 




CA 


-0.005(15)(20) 


-0.083(5) 


-0.013 




Cy 


-0.66(27)(02) 


-0.32(7) 


-0.028 




Ct 


+0.17(7)(1) 


N.A. 


-1-0.020 




bv 


+1.57(2)(1) 


M A 

iN.A. 


-l-l.lOo 




by^ 


+1 62f3)fl) 


1.54(2) 


-1-1.273 




bA ^ by 


-0.56(9)(8) 


N.A. 


-0.002 




bp - bs 


-0.09(9)(3) 


N.A. 


-0.066 




bp - Ia 


-0.20(9)(7) 


N.A. 


-f0.002 


bA 


^bp + bs/2 


+0.48(1)(2) 


N.A. 


-hO.584 


bA 


-bp + bs/2 


-H0.49(2)(2) 


N.A. 


-F0.673 




bA 


+1.01(09)(09) 


N.A. 


-M.104 




bp 


+0.81(14)(16) 


N.A. 


4-1.105 




bs 


+0.90(17)(13) 


N.A. 


-t-1.172 



Table 2 

Our best estimates for normalization and improvement 
coefficients, compared to previous results where available, 
and to tadpole improved 1-loop perturbation theory (us- 
ing uo — 0.8778). An estimate of the perturbative errors 
is ^ 0.02. The first error is statistical, the second an 
estimate of the O(a^) error. See text for details. 



a five-point discretization of 9^V^ when using the source 

j(21) ^^(21)_ 

We collect our best results in Table ^. These are ob- 
tained using the two-point derivative in SSi. The differ- 
ence between two- and three-point discretization is added 
as an additional error. Our results are compared to pre- 
vious non-perturbative results, and to those of tadpole- 
improved 1-loop perturbation theory. The latter are ob- 
tained using 1-loop results available in the literature ^ 
12 1 . We note that for the cx and the bx tadpole im- 
provement is equivalent to using the boosted coupling 

/uq in the 1-loop result. The same is not true for 



the Zx and bx- The conversion factor between the 
two definitions of the mass-dependent improvement co- 
efficients, bx/bx — Z\Zg/Zp, turns out to be very close 
to unity (see Table [|). Thus, with the exception of by and 
bA — bp + bs/2, which are determined quite accurately, we 
do not convert our results back to the standard definition 
bx- 
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We iraw the following conclusions. First, the method 856 (1998) . 



we have introduced appears to have practical utility. In 
the best channels, the statistical and systematic errors in 
the determination of the differences bA — by and bp — 65 
are small compared to the values, 6^ ~ 1, of the coeffi- 
cients themselves. Second, cy is determined rather poorly 
using our WI, and this accounts for a substantial frac- 
tion of the errors in Z^, Zp/Zg, and ct- For ca and cy, 
the Schrodinger functional method gives results with 
much smaller errors. Third, we have found, in some cases, 
substantial disagreement between the perturbative predic- 
tions and our non-perturbative results. The most striking 
cases are by and bA — by- These differences could be ef- 
fects of O(a^), and we intend to investigate this issue by 
repeating the calculations at different values of a. Finally, 
there are also differences between our results and those 
of the ALPHA Collaboration. We anticipate that these 
differences are due in part to our use of an action that is 
not fully 0{a) improved. We have initiated simulations at 
csw = 1-769 to verify this. 
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